Journal of Computational Neuroscience

We present a novel astro-neural-field population model with application to migraine-related cortical spreading depolarization. The model is composed of four spatio-temporal state variables: excitatory and inhibitory membrane potentials, astrocytic potassium uptake recruitment, and extracellular potassium concentration. Extending a previous neural field model, we incorporate activity-dependent astrocytic potassium clearance via a nonlinear term coupled to astrocyte dynamics. The astrocyte transfer function, like its neural counterpart, exhibits three regimes governed by extracellular potassium, capturing its effect on clearance. This yields a more comprehensive framework, better fits experimental data, and provides new insights into the mechanisms of cortical spreading depolarization.

Fly compound eyes pool signals from photoreceptors that sample the same region of visual space through neural superposition in the lamina. The optical axes of photoreceptors projecting to a single lamina cartridge are not perfectly parallel, but instead converge at a point a few millimeters in front of the eye. At short viewing distances (1-10 mm) this leads to distance-dependent differences in receptive field overlap. We explored whether it was possible that flies could sense depth in this "personal space" purely from the geometry of neural superposition. To this end, we combined a computational model of fly eye optics with a disparity-tuned lamina model originally developed for stereoscopic prey capture in praying mantises, and simulated the responses of lamina monopolar cells to moving stimuli at different distances. Across variations in stimulus parameters and lamina models, we found that lamina neuron responses indeed contain a distance-dependent component, which can overall be summarised as an enhanced response due to temporally overlapping receptor responses at a critical distance of 3-4 mm, the convergence distance of photoreceptor axes. Depending on the lamina model and stimulus size, either response amplitude or onset gradient, or both, exhibited this peak. We further show that changes in eye size systematically shift this preferred distance, such that larger flies had a peak response at a greater distance. Our results demonstrate that lamina cell responses may contain a robust, geometry-derived component that is specific to object distance and invariant to other stimulus properties. This suggests that neural superposition, beyond improving sensitivity, may function analogously to a light-field camera system that is effectively "focused" on a behaviorally relevant distance. Author summaryIn this study, we explored how neural superposition in fly compound eyes may lead to an enhanced response to objects at a behaviorally relevant distance. The optical axis of neurally pooled photoreceptors, from neighbouring ommatidia, are not parallel, but converge at a point a few millimeters in front of the eye, which should lead to the strongest correlation of their signals at that distance. This idea was tested by combining a geometric model of the fly eye optics with a computational model of how the lamina cells process the responses of photoreceptors, and evaluating the output for a moving bright dot at different distances. We found that the simulated lamina cell has its fastest response, as measured by the onset gradient, to a stimulus at the distance where the photoreceptors converge (3-4 mm). This distance scales with the simulated eye size and corresponds to behaviorally relevant distances for fly behavior. This suggests that neural superposition may act to enhance the response to objects at a critical distance in the early visual processing of flies.

This study presents improvements to the Hodgkin-Huxley (HH) models of ionic conductance and action potential generation. Sodium and potassium conductances are expressed by a single analytical formula describing the impulse response of a convolution of exponential distributions within a short-memory integration space. Treating transmembrane ion transit duration as a random variable, conductance profiles are interpreted as realizations of the probability density functions governing ionic movements. Applying the central limit theorem, the lognormal distribution emerges as the asymptotic profile of ionic conductances, constituting a fundamental primitive for such biosignals. A temporal state-transition paradigm describes the action potential waveform through four successive membrane potential transitions. Applied to electrophysiological recordings from lamprey reticulospinal neurons, this framework enables indirect estimation of key physiological quantities, including depolarization threshold, Nernst potentials, and net ion fluxes across the membrane. These advances open new perspectives for parameter estimation from experimental data and neuronal network simulation.

A novel mathematical framework to define the threshold of action potentials in excitable cells is presented. Unlike previously applied methods that rely on approximations or bifurcations, the approach focuses on the geometry of membrane potential trajectories. The changes in concavity during the upstroke of an action potential can be directly obtained from a time series of voltages. The concavity criterion is then extended to models based on autonomous dynamical systems where the changes in concavity can be obtained analytically from a curve of inflection points in phase space. The inflection point manifold defines a region required for excitability: all the orbits that cross it contain action potentials, and all the trajectories that contain action potentials are in it. This analytical principle can then be used to define excitability in a dynamical system, and also a measure of excitability that enables quantification and comparisons of excitability across dynamical system. The measure provides a way to compare the excitabilities of systems that model neurons with different electrophysiological phenotypes and consider different stimulus conditions. The traditionally vague physiological concept of electrical excitability is transformed into a rigorous analytical description by considering the time-dependent curvature of the membrane potential. The criterion is robust across smooth, single compartment models of electrical excitability and can be can be extended to single compartment models in higher dimensions, and multicompartment models as well.

Fixational eye movements (FEMs) comprise the involuntary small scale eye motion conducted during fixation on a stationary stimulus. As a consequence, the visual information can be spread across multiple photoreceptors reducing the local signal-to-noise ratio. Yet, the signals transmitted by individual photoreceptors adapt to constant stimulation so that an entirely still scene would eventually fade from view. Because FEMs convert a stationary stimulus in the world to a temporally varying one on the retina, they can act to prevent this stimulus fading. Thus, FEMs can be understood as a sampling protocol than needs to be adjusted to the underlying processing circuitry. We analyse the impact of FEMs on the rate of information acquisition at the level of the retina for two common tasks of the human eye that typically go hand in hand: detection and localization. Here, we build a simple analytical model of visual perception, i.e. we subject a continuous receptor array to a stimulus moving across the retina as a consequence of FEMs with receptor excitations depending on past stimulation through a linear response function. Using Bayesian inference we quantify both the probability of detection and the accuracy of localization as a function of parameters controlling eye movements and stimulus. We find that localization of a stimulus is equivalent to the detection of the stimulus gradient. This allows us to discern optimal properties of eye movements for the respective tasks and provides a link between two typical psychophysical observables: detection thresholds and Vernier acuity. Our analysis suggests that typical human FEMs tend to facilitate localization at the expense of detection. Simply put, if you can see a stimulus you also know where it is. Finally, we propose a variety of experimental protocols to investigate the interplay between FEMs, detection, and localization with the potential of inferring intrinsic properties of an individuals visual system.

Visual signals from the two eyes merge gradually as they pass through the primary visual cortex (V1). Here we use a computational model of Macaque V1 to study the first stage of this integration along the magnocellular pathway, in layer 4C, aiming to infer neuroanatomical origins of binocular response. It is known that neurons in layer 4C are predominantly monocular, though some do exhibit varying degrees of binocularity. We find (1) the emergence of narrow binocular strips along borders of ocular dominance columns (ODC), a finding that aligns with experiments; (2) most consistent with data is when 10 - 30% of interactions near ODC boundaries are cross-columnar; and (3) feedback from layer 6 is largely monocular. These results were obtained through systematic hypothesis testing using a multiscale model that is orders of magnitude faster than its biologically-detailed predecessors. We propose that multiscale modeling can be an effective tool for bridging anatomy and function.

Determining the amount of information a neuron receives per unit of time is key to understanding brain connectivity and how neural networks encode and transmit information. In particular, estimating this information flow by distinguishing between excitatory and inhibitory synaptic contributions is critical to understanding neural network function, as maintaining the excitation-inhibition (E/I) balance regulates neuronal excitability and circuit stability, whereas its disruption can lead to a plethora of brain disorders, including neurodegenerative and psychiatric conditions. However, because synaptic conductances cannot be measured directly, inverse methods are required to infer them from the membrane potential --a readily measurable quantity. Although partial solutions have been proposed, accurately estimating these conductances remains a significant challenge due to the complexity and diversity of the inputs. This is particularly true in the spiking regime, where neurons actively fire. In this work, we introduce a novel computational strategy that combines two critical metrics extracted from the time course of the membrane potential recording: the amplitude of the spike and the interspike interval. By using these quantities, the proposed method enables the accurate separation of excitatory and inhibitory contributions, yielding highly favorable results in the spiking regime. Author summaryQuantifying the continuous stream of inputs a neuron receives is key to understanding brain connectivity. Inside the brain, individual cells must maintain a tight balance between excitation and inhibition (E/I) to process information correctly, as any disruption in this equilibrium can impair its functionality. However, directly measuring the underlying excitatory and inhibitory synaptic conductances is technically challenging, and existing mathematical tools often fail when neurons enter their active firing regime. In this work, we introduce a novel computational strategy designed to extract and separate these time-varying conductances directly from the neurons spiking activity. By dynamically tracking just two accessible metrics - the amplitude of the spikes and the time intervals between them - our algorithm estimates both conductance profiles with high precision. Furthermore, we demonstrate that this procedure is highly robust against realistic experimental noise and data variability, providing an accessible framework that does not require complex hardware or an unfeasible number of repetitive experimental trials. By tracking changes in the E/I ratio of the synaptic input, this method provides an efficient approach to detecting pathological imbalances and understanding how local connectivity shapes cellular functionality.

Classical neurogeometric models describe the primary visual cortex as a fibered structure in which retinal position and local orientation are coupled through the geometry of the roto-translation group. We extend this approach to the visuomotor cortex by modeling it as an assemblage of visual and motor cortical geometries. The model combines orientation-selective representations, analogous to those of the primary visual cortex, with movement-direction-selective representations, analogous to those of the primary motor cortex, in order to describe the mixed visual and motor selectivity observed in the visuomotor cortex. We introduce a coupled visuomotor structure in which visual orientation and motor direction coexist over a common spatial plane and interact through a relative-orientation constraint. Neural responses are modeled by orientation- and direction-dependent profile functions, and preference maps are obtained from vectorized population responses. Numerical simulations generate visual, motor, and mixed visuomotor response maps. A competition rule between visual and motor responses produces incidence ratios close to experimental observations in macaque visuomotor cortex. This framework provides a first neurogeometric approximation of visuomotor functional architecture and a mathematical setting for studying visually guided action.

The awake thalamus is known to be able to filter primary sensory input with and without external modulation. Through the construction and analysis of a novel computational model of a lateral geniculate thalamocortical neuron, we demonstrate how the processing of sensory retinal input is influenced by the underlying thalamic dynamic state. Our model, using only currents verified against expression data from publicly available datasets, is the first to produce five experimentally established distinct dynamic firing regimes. We demonstrate that the thalamocortical cell transitions between these dynamic states in response to glutamatergic signals from the cortex or cholinergic arousal signals coming from the brainstem. We focus on signal processing in the model dynamic states associated with the awake thalamic alpha rhythm where we find that the ability of retinal inputs to generate thalamic spikes is a balance between the timing of retinal spikes, the excitability break imposed by the M-current, and the decay time of the L-type calcium current. Finally, we explore how these two currents help the thalamus process extra-retinal rhythmic inputs, showing the model produces entrainment to slower inhibitory and excitatory rhythms, as well as detailing the importance of nesting faster frequency rhythms within slow cycles for successful thalamic transmission. Our results suggest that the awake alpha rhythm is indirectly causal by acting as a marker for the interaction of these two currents. This biophysically-constrained lateral geniculate thalamocortical cell model generates predictions regarding rhythmic dynamics under different arousal states, thalamic control of retinogeniculate transmission, and the possible impacts neurological disorders, like schizophrenia, have on thalamic processing. Variations of this model could be used to explore the functions of higher order thalamic nuclei, thereby extending its use to investigating more complex cognitive processes. Author summaryThe thalamus generates multiple distinct brain rhythms, processes primary sensory inputs, and modulates its output using feedback signals. Previous computational models of the thalamus have typically focused on a subset of these three thalamic functions without drawing relationships among them. Here we present a novel computational thalamic cell model that unites these thalamic processes. We focus on the awake alpha rhythm, a well known thalamic oscillation, and show that it is a signature of a critical working state that enables the experimentally observed thalamic filtering of retinal signals. Additionally, we find this state is optimal for processing and passing non-sensory rhythmic signals. Our model generates testable predictions about which ionic currents control the transmission of external signals. It highlights the roles of two currents from our model that do not have specified functions in the awake thalamus in previous computational models. The work concludes with hypotheses about why neurological disorders that perturb the thalamus from this alpha rhythm working state lead to significant processing errors locally within the thalamus and globally within the brain.

Neural synchronization is fundamental to brain function and, when it becomes excessive, underlies pathological conditions such as epilepsy. Among brain regions, the temporal lobes, and the hippocampus in particular, exhibit the highest epileptogenic potential, with mesial temporal lobe epilepsy representing the most prevalent form of the condition in humans. Within the hippocampus, extracellular potassium dynamics are central to non-synaptic epileptiform activity, and astrocytic potassium buffering mechanisms have emerged as key regulators of network excitability. Yet the specific contributions of astrocytic gap-junction coupling and potassium spatial buffering to neuronal synchronization across different spatial scales remain poorly understood. To address this gap, we developed a microcircuit biophysical model consisting of two astrocyte-neuron modules, each comprising one astrocyte coupled to five neurons. Astrocyte-neuron interactions are mediated exclusively through shared extracellular potassium dynamics. Using a reduced astrocyte model that captures both local membrane and syncytial potassium buffering, we systematically investigated how astrocytic potassium handling shapes neuronal activity patterns and inter-module synchronization. Our results demonstrate that astrocytes prevent the emergence of pathological states -- such as sustained ictal activity and depolarization block, by stabilizing extracellular potassium levels. Furthermore, we show that astrocytic gap-junction coupling strength critically regulates phase synchronization between neuronal modules: stronger coupling promotes inter-module synchrony under physiological conditions, whereas impaired astrocytic function drives networks toward pathological hypersynchronization when extracellular potassium is elevated. These findings support the hypothesis that astrocytic networks impose modularity on hippocampal neuronal assemblies, and suggest that astrocytic connexins may represent a relevant therapeutic target in epilepsy and other disorders characterized by aberrant neural synchronization. Author summary

We have studied the emergence of slow relaxation oscillations in next generation neural mass models with spike frequency adaptation. Relaxation oscillations connect low firing state (Down state) to high firing state (Up state) via the slow adaptation. In the examined cases, the orbit relaxes towards the Up State via a sequence of collective damped oscillations (peaks of activity), thus revealing population bursting dynamics. The slower is the adaptation time scale the higher is the complexity (number of peaks) displayed by the relaxation oscillations. In particular, a chaos-induced spike-adding mechanism regulates the increase in the number of peaks. In analogy to what found in the Hidmarsh-Rose neuron model, two different types of chaotic behaviors have been identified: Population Spiking and Population Bursting Chaos. The increase of the adaptation strength leads to shorter (longer) Up (Down) state durations somehow mimicking the effect of charbachol in in vitro experiments, where spontaneous slow waves are observed. Indeed, the scenario depicted in [1], where an increase of the concentration of carbachol induces a transition from anesthesia-like to sleep-like dynamics is consistent with our results based on the variation of the adaptation strength. HighlightsO_LISpike Frequency Adaptation (SFA) promotes the emergence of Slow Relaxation Oscillations C_LIO_LISpike-adding mechanisms, controlled by SFA, lead to Relaxation Oscillations of increasing complexity C_LIO_LITwo types of chaotic behaviours: Population Spiking and Population Bursting Chaos C_LIO_LISFA regulates Up and Down States durations and their correlation C_LI

The activity of neural populations typically encodes more information about sensory or motor variables than can be captured by point estimates of the variables. We present and compare two approaches to quantifying this additional or ancillary information and its relationship to uncertainty: the mutual information between activity and estimation error, and the Fisher information loss, which can be interpreted in terms of curvature in information geometry. We show that deviations from Gaussianity of estimation errors, including the long tails frequently observed in human behavioural tasks, are an expected corollary of the presence of ancillary information. However, populations with similar distributions of estimation error can differ substantially in their ancillary information content depending on the noise characteristics. For a given population tuning and noise model, our results quantify an upper bound on the information about uncertainty that can be obtained from population activity alone: behaviour demonstrating knowledge in excess of this bound would indicate access to a separate source of information about uncertainty. Finally, we contrast the effects of external noise and decreasing internal signal strength on ancillary information and the Gaussianity of errors. Our work directly relates knowledge about uncertainty to non-Gaussianity in sensory estimates, and establishes a coherent theoretical foundation for investigating the basis of metacognition in neural population activity. Author summaryThe brain processes sensory evidence about the external world via inherently noisy neural activity. As a result, behavioural judgments - such as estimating the direction of a moving object - are fundamentally uncertain. While animals, including humans, routinely use uncertainty to guide decisions under risk, how neural populations represent this uncertainty remains unclear. In this work, we show how the same neural activity used to decode a sensory variable can also provide information about the estimates reliability. We introduce a mathematical framework to quantify this "ancillary information" directly from a neural populations encoding model. We demonstrate that ancillary information predicts non-Gaussianity in estimation errors and sets an upper bound on metacognitive sensitivity (how accurately subjective confidence tracks performance). Crucially, we show that neural populations with distinct noise characteristics can yield near-identical estimation errors while providing very different degrees of uncertainty information. This highlights the importance of evaluating ancillary information, not just error patterns, when comparing competing models of sensory coding.

Stabilizing synaptic plasticity together with the neurons activity has remained a central challenge in theoretical neuroscience since the introduction of Hebbian learning principles. Classical Hebbian learning rules typically lead to unbounded synaptic growth, motivating the development of stabilization mechanisms such as normalization methods, BCM-type learning, synaptic scaling and others. While these approaches can prevent divergence, they can also exhibit different limitations e.g. resulting in too-sparse synaptic configurations or leading to poor scalability with increasing network size. A recently introduced meta-plasticity mechanism, termed annealed linear learning (ALL), dynamically reduces the learning rate as neuronal output increases, thereby preserving stable and interpretable fixed-point behavior of the output. However, the original formulation leads to an irreversible decay of the learning rate, preventing adaptation to changing environmental conditions. To address this, in the present study, we balance learning rate reduction at large outputs with recovery at small outputs and in addition introduce forgetting that gradually reduces synaptic weights. These extensions allow the system to discard outdated representations and adapt to novel input conditions. Analytical investigations demonstrate that the favorable output fixed-point properties of the original ALL framework are preserved under the extended rule. Furthermore, simulations with an artificial agent show that the proposed mechanism enables robust and fast re-learning and adaptation in changing environments.

Natural physical environments constantly fluctuate across multiple timescales, often following a scale-free (1/f ) pattern where = 0.5 governs the fractional adaptation dynamics (Drew and Abbott 2006, Lundstrom et al. 2008). Here, we demonstrate how a multi-timescale sensory model successfully tracks these long-term trends to maintain stable encoding. Using an event-based Generalized Leaky Integrate-and-Fire (GLIF) paradigm, we found that a fast-adapting, single-exponential model with a short time constant{tau} [≤] 31.6 ms quickly crashes into complete refractory saturation when faced with large, low-frequency environmental shifts. In contrast, introducing a deep fractional memory tail of 1000.0 ms acts as an automated, high-pass balancing mechanism that continuously tracks and subtracts slow environmental variance. This predictive balancing prevents sensory collapse, anchors the mean firing rate to a steady homeostatic baseline, and maximizes coding efficiency for rapid, localized signals. Our results show that while a simple single-pole exponential model fails to retain history, a parallel bank of physiological relaxation processes converging on a target fractional profile t-0.5 provides the necessary historical memory to safely navigate natural stimulus fluctuations. Comfortingly, even a simplified three-pole approximation captures the bulk of this homeostatic benefit, making efficient fractional adaptation biologically viable at the sensory periphery without requiring infinite historical storage.

Biological neurons can perform nonlinear computations within their dendrites and support branch-localized plasticity. This raises the possibility that single cells can store memories more efficiently and with less interference by confining synaptic modifications to specific dendrites. We study a parallel-dendrite model performing online familiarity detection and compare three dendrite-update rules during learning: (i) independent thresholding, (ii) an interacting rule that adapts the target local dendritic activation per item, and (iii) an interacting n-winners-take-all (WTA) rule that constrains the number of updated branches per item. The interacting rules substantially improve capacity by limiting variance in memory responses and decorrelating weights across branches -- even when inputs are strongly correlated. These results suggest that competition among dendrites, consistent with resource-limited plasticity mechanisms, can enhance single-cell memory beyond non-interacting schemes.

Cortical populations exhibit a wide range of tuning properties, raising the question of whether such variability is a feature or a bug of cortical function. Prior work has shown that tuning diversity can improve population codes by mitigating the effects of correlated noise and increasing the discrimination and identification capacity of geometric representations. Motivated by these findings, we study a model in which a heterogeneous family of tuning curves, coding for a circular variable, is replicated at equally spaced preferred angles. We show that this heterogeneous population achieves better discrimination and detection than an equally sized homogeneous population constructed from shifted copies of the family's mean tuning curve, while using the same spike budget. Thus, homogeneous tuning is unstable under perturbations that preserve the mean tuning curve, because such perturbations leave metabolic cost unchanged while improving coding performance. We propose that such instability creates evolutionary pressure toward heterogeneity of tuning, making its prevalence a consequence of a process that optimizes coding performance under metabolic constraints.

Segregation of speech into syllables is a key step in neural speech processing. It relies on the alignment of neural activity with the rhythmic structure of speech. Two competing hypotheses explain this neural speech tracking, phase-resetting and evoked responses. While phenomenological modelling of these hypotheses has been successful, we still lack understanding of the underlying cortical circuits. To investigate these mechanisms, we evaluate whether a biophysical next-generation neural mass model can reproduce several features of neural speech tracking, using phenomenological models of the competing hypotheses as algorithmic baselines. We investigate the models dynamics with four tests: recreating in-silico an EEG experiment that identified a correlation between tracking strength and phoneme sharpness, computing the Phase Concentration Metric, testing the effect of varying syllabic rates, and evaluating the Inter Event Phase Coherence across phoneme onsets. While all of the models that we study reproduce the sharpness-tuned rhythmic speech tracking, the evoked model requires a pre-processed acoustic edge impulse stimulus. We demonstrate that the neural mass model is performing thresholded phase-resetting triggered by sharp onsets in the continuous speech envelope. This produces cross-frequency nested oscillations that qualitatively match an experimentally-observed dual-peak signature in the Inter Event Phase Coherence. Our results indicate that the biophysical neural mass model provides a mechanistic bridge between generic oscillatory dynamics in cortical populations and the cognitive computations of speech tracking. Indeed, the non-linear dynamics of the neural mass model offer an explanation for how peak-rate event representations in auditory cortex activity arise in response to continuous acoustic input. Significance StatementSyllable segregation is crucial but challenging as natural speech lacks clear boundaries, yet humans perform this computation effortlessly. Speech aligns neural activity to syllabic rhythms, predicting syllable timing, but the underlying cortical mechanisms remain unknown. Relating this macroscopic behaviour to neurobiology is challenging; however, next-generation neural mass models promise to resolve this. We demonstrate that these models reproduce sharpness-tuned tracking and acoustic edge extraction. Dynamical analyses indicate this occurs through thresholded phase-resetting to phoneme onsets, triggering cross-frequency nested oscillations. Our results both advance biophysical understanding of syllable segregation and validate the models capacity for simulating macroscopic neural activity. These models offer a bridge between the neurobiology of the auditory cortex and speech processing dynamics that phenomenological models cannot provide.

Deep brain stimulation (DBS) has been demonstrated to be a successful therapeutic intervention for neurological disorders, yet the mechanisms underlying its effects on neuronal circuits remain incompletely understood. In this study, we propose a comprehensive phenomenological computational model that accounts for the impact of electrical stimulation parameters on neuronal circuits while incorporating experimentally-validated synaptic and cellular constraints. We investigate how DBS pulses modulate spiking activity in populations of homogeneous neurons representing stimulated nuclei, systematically examining the influence of circuitry architecture, including synaptic connectivity strength (weak vs. strong) and organization (sparse vs. rich). To characterize how DBS-modulated neuronal activity propagates through downstream networks, we develop a simple encoder that reveals distinct encoding patterns arising from different architectural configurations of stimulated nuclei. Furthermore, by connecting stimulated nuclei to recurrently connected neuronal populations, we examine the propagation of DBS-modulated neuronal synchrony across various circuit motifs. Our results demonstrate that three critical factors shape DBS-modulated neuronal activity: (a) the intrinsic synaptic and cellular properties of stimulated nuclei, (b) the architectural organization of stimulated nuclei in terms of synaptic strength and connectivity density, and (c) the circuit motifs formed by postsynaptic targets of stimulated nuclei. This unified model provides a mechanistic framework for understanding DBS representation and propagation in neuronal networks, offering insights that may inform optimization of stimulation parameters for clinical applications.

The synaptic tagging and capture (STC) hypothesis explains how early-phase plasticity is converted into its late phase through the coincidence of synaptic tagging and plasticity-related protein (PRP) availability. Yet the biophysical basis of this process remains poorly understood. Based on the hypothesis that the interaction of actin and spine geometry implement the synaptic tag, we here investigate the associated PRP capture mechanism. We propose that capture is implemented by PSD remodelling which is gated by local membrane curvature at the PSD periphery. Using computational modelling, we show that curvature variations around the PSD that arise from long-term potentiation (LTP) inducing stimuli indeed enable a PSD growth, reproducing late-phase potentiation and the maintenance of structural LTP. We further explore how the timing of PRP availability relative to tag formation and the initial spine size determine the extent of PSD enlargement, yielding outcomes consistent with experimental findings. Hence, our results support a structural interpretation of synaptic tagging and capture in which a transient, actin-driven geometric state of the spine encodes the tag, and curvature-mediated PRP recruitment stabilises synaptic changes, and thus render spine geometry as a key biophysical regulator of memory consolidation.

A complete understanding of population coding requires connecting multiple levels of neural processing: individual responses, population representations, and behavior. We link these by relating the distribution of neuronal tuning properties to a populations representational geometry and its efficiency for perceptual tasks. We use theory, analysis of recordings from macaque primary visual cortex (V1), and simulations to reveal how diversity of tuning amplitude and bandwidth enhances the population code for visual discrimination and identification. Both types of diversity drive different, but complementary changes to the representational geometry. Amplitude diversity increases the Euclidean distance between the responses to different stimuli, while bandwidth diversity creates a larger angular distance between them. The first utilizes the range of firing rates available to neurons, and the second exploits the high-dimensional nature of population responses. Population codes can be improved using these two different geometric changes, and amplitude and bandwidth diversity provide biological mechanisms for doing so. HighlightsO_LI- Perceptual performance is improved both by increased diversity of response amplitude and increased diversity of tuning bandwidth. C_LIO_LI- Both kinds of diversity improve visual discrimination and identification. C_LIO_LI- Amplitude diversity improves discrimination more, and bandwidth diversity improves identification more. C_LIO_LI- Representational geometry reveals the mechanisms of these effects. C_LI